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Abstract 

It has previously been pointed out that the coexistence of infinite-range and short-range inter- 
actions causes a system to have a phase transition of the mean- field universality class, in which 
the cluster size is finite even at the critical point. In the present paper, we study this property in 
a model of bistable molecules, whose size changes depending on the bistable states. The molecules 
can move in space, interacting via an elastic interaction. It is known that due to the different 
sizes, an effective long-range interaction between the spins appears, and thus this model has a 
mean-field type of phase transition. It is found that the scaling properties of the shift of the criti- 
cal temperature from the pure short-range limit in the model with infinite-range and short-range 
interactions hold also in the present model, regarding the ratio of the size of the two states as a 
control parameter for the strength of the long-range interaction. By studying the structure factor, 
it is shown that the dependence of the cluster size at the critical temperature also shows the same 
scaling properties as a previously studied model with both infinite-range and short-range interac- 
tions. We therefore conclude that these scaling relations hold universally in hybrid models with 
both short-range and weak long-range interactions. 
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I. INTRODUCTION 



The divergence of the correlation length at the critical point is considered to be one of 
the most important properties of second-order phase transitions. However, it is also known 
that in phase transitions of infinite-range interacting systems that belong to the mean-field 
universality class, the correlation length does not diverge, and the spatial configuration is 
uniform with no domain structures or clustering. In a previous work, we studied a hybrid 
model with both short-range and weak long-range interactions.- The Hamiltonian of that 
model is given by 

U = (1 - a)Hjs + oHmf, (0<a<l) (1) 

with 

%s = -J^Wj, ( 2 ) 

(hi) 

and 

^mf = -^J2 a i a i = -^fX^r ( 3 ) 

i,3 hi 

We defined the model on the square lattice with periodic boundary conditions. Here, the 
strength of the infinite-range interaction is controlled by varying a. When a = 0, the system 
is equivalent to the pure short-range Ising model, and the system with a = 1 is equivalent to 
the pure infinite-range interaction model. The critical temperature in the previous hybrid 
model ([I]) shows a crossover from that of the pure short-range Ising model to that of the 
infinite-range interaction model. It should be noted that even for infinitesimally small a, 
the phase transition belongs to the mean-field universality class, and at the critical point, 
the spin configuration is uniform with no large-scale clustering. A scaling formula for the 
a-dependence of the critical temperature T c is found, such that 

T TmJIt^ ~ L773517c ^ = 1-773517^. (4) 

c c 

The correlation length £ c at T c is also found as a function of a, 

U<x,L) = Lf(La$) = Lf(La$) (5) 

for small a and large L. The function f(x) is a scaling function which asymptotically 
approaches \ jx for large x. 
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Although the long-range interaction in this model is rather artificial, recently it has been 
pointed out that spin-crossover materials and related materials show a similar kind of long- 
range correlation.- Spin-crossover materials are molecular crystals, in which the molecules 
can exist in two different states: the high-spin (HS) state and the low-spin (LS) state. The 
HS state is preferable at high temperatures because of its high degeneracy, while the LS state 
is preferable at low temperatures because of its low enthalpy. In addition to temperature, 
pressure changes and light exposure also often induce a phase transition in spin-crossover 
materials. Spin-crossover and related materials are used in many applications, because of 
their inherent bistability that leads to changes in optical and magnetic properties, etc.— 
Phase transitions in spin- crossover and related materials have been studied extensively in 
chemistry^— and recently also in physics.— ~— For a wide variety of applications, it is of 
great interest to study the ordering process in spin-crossover materials. 

This type of materials are also regarded as fundamental models for inter-molecular short- 
range and elastic interactions due to the lattice distortion, and an elastic interaction model 
has been proposed.- 1 ^ An important characteristic of this model is an effective long-range 
interaction due to the lattice distortion caused by the size difference between the HS (large) 
and LS (small) molecules. Even at the critical point, there exists no large-scale domain 
structures. However, in the elastic interaction model, there does not exist inter-molecular 
short-range interaction, so the model does not show any crossover from an effective short- 
range interacting system to an effective long-range interacting system. Therefore, developing 
quantitative hybrid models with both short-range and the elastic interactions for the critical 
behavior of such materials is very important. 

In this paper, we focus on a system with both elastic interactions and short-range Ising 
interactions. We perform Monte Carlo (MC) simulations to study the properties of this 
model at the critical point. The rest of the paper is organized as follows. In Sec. [Til we 
introduce a new model with both elastic and short-range Ising interactions, and we propose 
a relation between this model and the model studied previously- We find that the scaling 
formulae for the critical temperature and correlation length obtained in the previous paper- 
apply to this model, as well. In Sec. we briefly review the MC algorithms and give the 
result of MC simulations for the critical temperature, confirming the relation. In Sec. IIV| 
we similarly confirm the scaling form of the correlation length at the critical point. In Sec. 
IV] we summarize our results, and in appendix [A] we discuss in detail how we calculate the 
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correlation length at the critical point. 



II. MODEL 

In this paper, we adopt the following model with both lattice distortion and inter- 
molecular short-range interactions on the square lattice with periodic boundary conditions. 
From now on, L denotes the number of molecules along an edge of the lattice, so the total 
number of sites is N = L 2 . The Hamiltonian is 

H = "His + "Hnn + "%nnn + Heft, (6) 



with 



%IS — —JJ2{i,j) 

n nn = f E (i>i > [h -rj\- {Riin) + RMM 2 > 

Hnnn = ¥ E((L m )) |»*/ ~ T m \ - \/S (Rl((Tl) + R m ((T m )) 



(7) 



U e{i = (p-bflogg) Ef^i, 
where r« represents the continuous coordinate of the molecule i, and Ri(<Ti) is the radius of 
the molecule i which depends on the molecular state <Tj. Here <7j = +1 and —1 represents 
HS and LS state, respectively. Hereafter, we simply call Cj the spin state. His is the 
short-range pure Ising model. "H nn and "H nnn denote the elastic interaction Hamiltonians of 
nearest-neighbor and next-nearest-neighbor ((l,m)) pairs, respectively, and k\ and hi 
are the corresponding spring constants. The next-nearest-neighbor interaction is introduced 
to maintain the shape of the lattice, and the strength of k-i is not important as long as the 
global shape of the square lattice is kept. Here we take k 2 = fci/10. We define the pure 
elastic interaction model as 

"^Elastic = ^nn + ^nnn- (8) 

The molecular radius is determined by the local spin state : i?n for the HS state (large, 
(Tj = +1) and i?L f° r the LS state (small, (7j = —1). When |r~ Vj\ is equal to the sum of the 
radii Ri(<Ji) +Rj(aj), the corresponding contribution to the elastic energy has its minimum. 
In H e ft, which represents the ligand field, D denotes the energy difference between the HS 
state and LS state, and g denotes the ratio of the degeneracies of the HS state and LS state. 
The order parameter of these models (jHJ) and (JHJ) is defined as 

m ss ^p-, (9) 



which is related to the fraction of HS molecules, /hs, as m = 2/hs ~ 1- 

In order to see the competition between the short-range interaction and the lattice dis- 
tortion due to the molecular size difference, we consider the model ([H]) along the coexis- 
tent line, (to) = 0, as we studied in the previous work for the Ising model. Thus we set 
D — logg = 0. For simplicity, this situation is described by the present model §§§ in 

which D = and g = 1. In this model, a ferromagnetic second-order phase transition takes 
place at the critical temperature T c . Below T c , there exist two different ordered states, the 
HS state and the LS state. 

It has been found that in the pure elastic interaction model (jSJ), the spin configuration is 
uniform at the critical point, with no large-scale clustering.- The phase transition belongs 
to the mean-field (MF) universality class, and the spin correlation function approaches a 
non-zero constant in the long-distance limit. 



FIG. 1: (Color online) Typical configuration of the mixed elastic and short-range interaction model 
at the critical point. R^/Ri, = 1.1, for which T c = 0.54575. Light gray (red online) filled circles 
denote HS state spins, and dark gray (blue online) filled circles denote LS state spins. 

The origin of the long-range interaction in the elastic model is lattice distortion due to the 
size difference between the HS and LS states. We expect that the ratio of the radii R^/Ru 
controls the strength of the effective long-range interaction. Namely, R^/Rn plays a role 
similar to a in the previous hybrid model (pQ). The strength of the long-range interaction 
is given by the elastic energy, which is of the order of (Jfei/2)(-R H - #l) 2 oc (1 - f^) 2 . We 
therefore consider that the parameter to indicate the strength of the long-range interaction 
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is given by 



Rr 2 



a oc ( 1 - ) . (10) 



With this relation, we expect that the relations (j4]) and (j5j take the following forms, 




(A^,L)=Lf\Hl-^) ], ,12; 



and 

- Rw 

respectively. As in our previous work, 1 as long as the value of 1 — R^/R^ is small, we expect 
these formulae to be correct. In Fig. HJ we depict a typical configuration of the elastic and 
short-range interaction model at its critical temperature. Unlike the pure short-range Ising 
model, the spin configuration is uniform with no large-scale clustering, even at the critical 
temperature. A lattice distortion also occurs, and we observe an uneven system surface. 



III. CRITICAL TEMPERATURE 



In this section, we perform MC simulations to test the scaling relation for the critical 
temperature (TTTT) . For the simulation, we adopt the NPT-MC method^ for the isothermal- 
isobaric ensemble with the number of molecules N, the pressure of the system P, and the 
temperature T. In this paper, in order to exclude other effects than those due to the elastic 
interaction through distortion, we fix P = 0. We also fix the spring constants as k\ = 40 
and &2 = 4 as in our previous work.- The critical temperature of the pure elastic interaction 
model ([8]) is T c Elastlc ~ 0.2 for i?n — 1-1,- and we choose J = 0.1 in order to keep the two terms 
in the present model ([6]) of comparable magnitude near the critical temperature. Therefore, 
T C IS = 0.2269 • • • in these units on the square lattice. We use a standard Metropolis method, 
adopting periodic boundary conditions. In most cases, we performed eight independent runs 
of 4,000,000 Monte Carlo steps per spin (MCSS) for the each data with 100,000 MCSS for 
the initial equilibration. We confirm that the statistical errors are smaller than the marks in 
the following graphs. We fix R L = 1.0 and choose R u = 1.005, 1.008, 1.010, 1.015, 1.02, 1.05, 
and 1.1. 

We have previously pointed out that even infinitesimally weak long-range interactions be- 
come dominant in the thermodynamic limit.- In the case of the previous hybrid model (CD), 
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this property is explained by the fact that in a well coarse-grained Hamiltonian, the long- 
range interactions become stronger than the short-range ones. We need systems sufficiently 
large that clusters caused by the short-range interactions can be regarded as block spins.- 
Here we assume that this size dependent crossover phenomenon also takes place in the 
present model fl§]). 

Here, in order to estimate the critical temperature, we adopt a method we also used in our 
previous paper. We use the crossing point of the forth-order Binder cumulant^ U^L) for 
different system sizes to estimate the critical temperature T C {R^/ Rh) with high precision. 
The cumulant is defined as 

ML) s 1 - g|, (13) 

where m = 1/ L 2 Y,i^i- In the case of the pure elastic interaction model (JSJ) on the square 
lattice with periodic boundary conditions, the fixed-point value of the cumulant 2 is the same 
as the exact value for the infinite-range interaction model (j3J), ~ 0.27 ■ • •.— ^ With the 
radii ratio R^/Rb. = 1, the present model (jSJ) is equivalent to the pure Ising model on the 
square lattice with periodic boundary conditions. For this case the fixed-point value of the 
cumulant is the same as the value for the Ising model (j2J), ~ 0.61 • • • on the square 
lattice.— Other shapes of the system, boundary conditions, and anisotropy may lead to 
different values of U4, at the crossing point.— "22 



S 




0.4 

0.3 
0. 



o 


L=10 


□ 


L=20 




L=40 


V 


L=80 





L=160 



250 



0.255 




H 
t3 



0.6 



0.5 



0.4 



0.3 




o 


L=10 


□ 


L=20 




L=40 


V 


L=80 



0.265 



0.270 



0.33 



0.34 



0.35 



(c) (d) 

FIG. 2: (Color online) Temperature dependence of the Binder cumulant Ui{L) for (a) i?H = 1-005, 
(b) i?H = 1.01, (c) Rh = 1.02, and (d) Rn = 1.05. Points denote Monte Carlo data, and the solid 
lines are polynomial fits. The upper and lower horizontal lines denote the fixed-point values for the 
Ising model (Ul s ~ 0.61) and the infinite-range model ~ 0.27), respectively. The left vertical 

line in (a) represents the critical temperature of the pure Ising model. The linear system size L 
is 10, 20, 40, 80, 160, and 320 for circles, squares, upward triangles, downward triangles, diamonds, 
and asterisks, respectively. 



In Fig. [21 we see that the crossing points of the Binder cumulant decrease toward the 
mean-field fixed-point value, ~ 0.27- • — ^ from the Ising fixed-point value, Ul s ~ 

0.61 • • — as L increases. This indicates that a size dependent crossover occurs, and that 
the critical point of this model belongs to the mean-field universality class. We estimate the 
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FIG. 3: (Color online) Size dependence of the crossing point of the Binder cumulant (a) and an 
enlarged view of the same for small i?n (b). From above to below, we show simulation results with 
L = 20,40, and 80 for R n = 1.10 and 1.05, L = 20,40,80, and 160 for R n = 1.02,1.015,1.01 and 
1.008, L = 20,40,80, 160 and 320 for R H = 1.005. The points at l/L = are linear extrapolations 
from the two smallest nonzero values of l/L. Filled squares, filled circles, diamonds, downward 
triangles, upward triangles, open squares, and open circles represent T C (L) for i?n = 1-10, Ru = 
1.05, i? H = 1.02, R u = 1.015, i? H = 1-010, i? H = 1.008, and R L = 1.005, respectively. 

critical temperature as follows. Assuming that the critical behavior of the model belongs to 
the mean-field universality class, we get a series of upper bounds on the critical temperature 
as the temperature at which U^{L) crosses L 7 | 1F . Lower bounds are given by the cumulant- 
crossing temperatures of U^{L) and U^{L/2). In Fig. [31 bars denote those upper bounds and 
lower bounds of the critical temperature, and we plot the middle points of those bounds 
by bullets. Increasing L, the temperature range between the upper bounds and the lower 
bounds becomes narrow. For each value of R^/Ru, we extrapolate those middle points for 
the two largest system sizes versus l/L to obtain the corresponding critical temperature in 
the thermodynamic limit. With this method, we expect that a critical temperature of a 
hybrid system with both short-range and long-range interactions can be obtained accurately 
in general. The data collapse well onto a straight line as shown in Fig. HI confirming the 
scaling form ( fTTl) for i?H smaller than 1.02. The points fall above the line for large values of 
1 — R-^/Rn. In this model, we fix the strength of the Ising interactions, while in the previous 
hybrid model ([1]), the factor 1 — a multiplying the Ising Hamiltonian causes a deviation of 
the critical temperature in the opposite direction (see Fig. 4 of Ref. 1). 
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FIG. 4: (Color online) The ( 1 — J dependence of the shift of the critical temperature in a log-log 



plot. The oblique dashed line (blue online) represents a numerical fit of the 4 leftmost points to 



IV. CORRELATION LENGTH 

In this section, we perform MC simulations to test the scaling relation for the correlation 
length at the critical point ( jl~2l) . We previously calculated the critical correlation length 
in the previous hybrid model ([I]) from MC simulations.- In that paper, we used the spin- 
correlation function c(r) = (a(r')cr(r' + r)) to estimate the correlation length, excluding the 
contribution of long-range correlations from the correlation function. In the present model 
(EJ), because of strong anisotropy of the correlations,- it is not practical to calculate the 
correlation length £ (Rl/Ru) with this method. 

Here we instead obtain the correlation length from the structure factor, 



which is readily measured in scattering experiments. All the contributions of the long-range 
correlations are given by S(k = 0), so we can easily exclude them from the calculation. 
Here we note that there exist many experimental studies on the structure factors of mag- 
netic materials, 2 - obtained by neutron scattering. In the case of spin-crossover materials, 
the structures have been studied by single-crystal x-ray diffraction experiments.— 1 ^ The 
structure factor for the HS/LS state domains discussed in the present paper can be obtained 





l,m 



(14) 
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from diffuse x-ray scattering. 

In the pure short-range Ising model at the critical temperature, the structure factor has 
its peak at k = with infinitesimally narrow width in the thermodynamic limit. However, 
in the hybrid model with both short-range and long-range interactions, the system prefers 
a spatially uniform configuration and the cluster size is suppressed even at the critical 
temperature, and the peak has a finite width in the thermodynamic limit. Here we obtain 
the correlation length at the critical point, £ c by calculating the characteristic peak width 
of the structure factor. 
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FIG. 5: (Color online) The structure factor of the pure short-range Ising model (circles) and the 
hybrid model with i?n = 1-02 (squares), at each its critical point, (a) L = 160 and (b) L = 20. 



In Fig. El we depict the structure factor of the pure Ising model at its critical temperature 
(circles) together with that of the hybrid model at its critical temperature (squares). Figure^] 
(a) shows the structure factor of the system for L = 160, and Fig. [5] (b) is for L = 20. We 
find a qualitative difference between these figures. Namely, we find a flat region in S(k) of 
the hybrid model around k = in Fig. [5] (a), while S(k) for both the pure short-range Ising 
model and the hybrid model show similar peaks in Fig. [5] (b). 
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The structure factor of the hybrid model consists of two parts: A 5-function at k = 
and a broad peak reflecting short-range order with the correlation length £ c . The two-parts 
are superimposed in Fig. (a). The flat region belongs to the diffuse peak, which is due to 
the finite cluster structure of the hybrid model. We find that for sufficiently large systems, 
the diffuse peak and the 5-peak at k = are well distinguished. In those systems, we define 
the characteristic peak width |fc pea k| as the spectral peak width of the diffuse peak as shown 
by the arrows in Fig. [5] (a). 

In Fig. [5] (b), for small L, the resolution in fc-space (2tt/L) is rather coarse, so the data 
points do not reflect the finite width of the peak in the hybrid model well. It is hard to 
distinguish between the structure factors of the Ising model and the hybrid model in such 
small systems. 

We note that the sum rule of the structure factor is J2k S(k) = N. In the hybrid model, 
S(k) for large \k\ is larger than that of the pure short-range Ising model because S(k) 
around k = is suppressed in the hybrid model. 

We can estimate the correlation length £ from the structure factor by considering the 
first moment of k' 1 , 

r 1 )^^), as) 

where k = \k\. In our present model, the correlation function of the spin configuration is 
given by2>22 

= c(r M ) = c SR (\r t - r,|) + c LR (16) 

Here, c SR (\ri — rj\) denotes the contributions from the short-range interactions, and c LR 
denotes those from the long-range interactions. We note that the long-range correlation c LR 
does not depend on the distance |r*j — r 3 -|. In Eq. (fT5l) . the contribution from the long-range 
correlations c LR (i.e., S(k = 0)) is excluded, so we use the asymptotic formula for the pure 
short-range Ising model,— 

C M « -^r^. (i7) 

In the thermodynamic limit, we replace the sum in ( |T5|) by an integral: 

poo pin 1 poo p2ir 

(AT 1 ) = / / kdkdcf)- / / rdrd6e ikrcose c(r) (18) 

Jk=2n/L J(f>=0 k Jr=0 J 8=0 

/•oo roo 

/ dk / drr^e^^Mkr) (19) 

Jk=0 Jr=0 

/•oo poo 

^-' n / / R- 71 e~ R J (K)dKdR (20) 
Jk=o Jr=o 
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oc e~ v = e 3/4 - 



(21) 



In the last line, we use the substitutions, K = kr, R = r/£, and Jq{K) is the Bessel function. 
There is no singularity in the integral, and (A; -1 ) gives a power of the correlation length. 
Therefore we define the critical correlation length £ c at the critical temperature T c as follows, 



where (k 1 ) c represents the value at the critical point. The scaling relation for (k 1 ) c is, by 
making use of (JT2"|) . 



where g(x) oc - for x — > oo. 

We obtained (A; _1 ) c and the corresponding values of £ c for various values of L and R^/R^. 
Those are plotted in Fig. [61 In the enlarged view of Fig. [6] for small L, we find inflection 
points. For small L, £ c (Rl/Rh, L) grows faster in the hybrid model than in the pure Ising 
model. 

Here let us consider the reason for this behavior. This fast growth is due to the coarse 
resolution in fc-space. The characteristic peak width of the hybrid model does not represent 
well the structure factor in the small systems (Fig. [5(b)), and the correlation length of the 
hybrid model is overestimated in Eq. (|23|) . This extra size dependence for small size comes 
from the discontinuous nature of fc-space, and is not essential for the present purpose. 

For large L, the growth of £, c {Rl/ Rn, L) becomes much slower than in the pure Ising 
model. For i?n = 1-02 and L = 320, ^(R^/Ru, L) is almost saturated. For other values 
of i?H, £c(-Rl/-Rh, L) is expected to saturate for still larger systems. It means that for 
sufficiently large systems, the difference of the structure factor of the hybrid model and the 
pure short-range model is clear. 
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(22) 




(23) 
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FIG. 6: (Color online) The size dependence of (A; _1 ) c 3 oc £ c at the critical point T^R^/Rh) (a), and 
an enlarged view of the same for small L (b). L = 2, 4, 6, 8, 10, 14, 20, 30, 40, 80, 160 and 320. Circles, 
squares, upward triangles, downward triangles, diamonds and crosses represent T = 0.23107 for 
R u = 1.003, T = 0.23428 for R u = 1.005, T = 0.23965 for i? H = 1.008, T = 0.24345 for i? H = 1.01, 
T = 0.2537 for i? H = 1-015 and T = 0.2649 for R L = 1.02, respectively. The dashed line (blue 
online) represents T = 0.2269 for Ru = 1.00, the pure Ising-model limit. 

The same behavior is found in the previous hybrid model ([T]). In both models, the 
previous hybrid model (JTJ) and the present model ([6]), the long-range interactions suppress the 
clustering of spins. Thus there appear tightly correlated effective block spins. In the coarse- 
grained Hamiltonian with the block spins, the long-range interactions become effectively 
stronger than the short-range interactions.- 

According to this picture, for L smaller than the clusters caused by the short-range 
interaction, the long-range interaction is irrelevant. Thus, for small L, the system is an 
effective short-range interacting system. On the other hand, for sufficiently large L, the 
system is an effective long-range interacting system. We refer the crossover length as L c \, 
which we define as L c \ = l/|fc pea k|- The inflection points in Fig. E] also indicate the crossover 
length of the system size, L c \. For L > L c \, the correlation length at the critical point, 
£ c (-Rl/-Rh, L), tends to be saturated. For L < L c i, £ c (-Rl/-Rh, L) grows faster than linear. 

We also found that the peak position of x(L) = j^((M 2 ) — (|M|) 2 ) shows similar behavior. 
The peak position, the effective 'critical point' for a system of size L , saturates at the critical 
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temperature in the thermodynamic limit: 



> -ixh / V -tin 



(24) 



As in the previous work,- we found a non-monotonic dependence of the peak position as 
a function of L. For small L, the peak positions approach T c from the high-temperature- 
side, while for large L, the peak position is on the low-temperature-side and eventually 
approaches T c from below. This non-monotonic behavior indicates the crossover from the 
effective short-range system to the effective long-range system. It is an interesting problem 
to study the relation between the crossover phenomena in x(L) and the correlation length 
£ c in Fig. [HI We expect that this general system size dependence occurs in any hybrid model 
where effective spins play an important role. 

In appendix A, we discuss in detail how the correlation lengths are measured. In Fig [71 
we plot the data for L > L c i(_Rl/-Rh) in a scaling plot of the form (1231) . We find that the 
data collapse onto a scaling function and thus we conclude that (123|) is justified. 

^ 1 
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FIG. 7: (Color online) Scaling plot of the correlation length at the critical point. Circles, squares, 
upward triangles, downward triangles, diamonds, and crosses represent T = 0.23107 for i?n = 
1.003, T = 0.23428 for i? H = 1-005, T = 0.23965 for R n = 1.008, T = 0.24345 for R u = 1.01, 
T = 0.2537 for = 1.015, and T = 0.2649 for i?L = 1-02, respectively. Data are included only 
for L > L c \. The dashed line is proportional to y = The data are in good agreement with 

the scaling relation, but they converge more gradually to the asymptotic formula than when the 
correlation length is calculated from the correlation function. 1 
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V. SUMMARY AND FURTHER DISCUSSION 



In general, spin-crossover and related materials have both short-range and long-range 
interactions. If the elastic potential is chosen spin-state dependent^ 1 ^, in which different 
potential functions (coupling constants) are given for LS-LS, LS-HS, and HS-HS molecular 
pairs, it is considered that the interaction has intrinsically both short-range and long-range 
components. In this paper, we have studied critical properties of the present model ([6]) with 
elastic and short-range interactions. By Monte Carlo simulations we confirmed that the 
present model obeys a scaling relation for the shift of the critical temperature as a function 
of the strength of the long-range interaction, which is determined by the ratio of the radii 
of the molecules in both states. We similarly confirmed that the present model also obeys 
a scaling relation for the correlation length at the critical point, which was found in our 
previous work for the hybrid model with infinite-range and nearest-neighbor interactions.- 1 
Although the origin of the long-range interaction is very different in the previous model and 
the present model, we found several kinds of universality in their critical properties. 

In the present model, because of the anisotropy in the correlation function, we estimated 
the correlation length from the structure factor. By this method, the correlation length can 
be measured in scattering experiments on real materials. Because spin- crossover materials 
usually undergo first-order phase transitions, experiments on the structure factors of these 
materials have been so far only done for HS/LS ordered phase.— ^ Our results show that 
the structure factor of spin-crossover materials should exhibit a peak with finite width at the 
critical point due to the elastic interaction. We hope such unique behavior will be observed 
by single crystal x-ray diffraction along the coexistence line and at the critical point. 

It was confirmed that the present model ([6]) possesses an effective long-range interac- 
tion. The details of the effective long-range interactions introduced by the elastic degrees 
of freedom in the present model are not known. Only for d — 1,— it has been shown rigor- 
ously that the model can be mapped onto an short-range Ising model. There has also been 
much previous research on three-dimensional elastic solids, and it is generally argued that 
the dominant long-range interactions are of a dipole-dipole nature, ~ 1/r 3 . 26 Although the 
effects of distortions in the present model are not identical to those in the classical elastic 
media for which these results were obtained, we assume that the elastically mediated inter- 
actions in our model also are of such long-range type. However, in our present model the 
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spin correlation function shows infinite-range correlations above the critical temperature.- 1 ^ 
Therefore, we assume that there are similar but unknown infinite-range interactions in the 
present model. Understanding the mechanism by which the infinite-range interactions arise 
remains an intriguing problem for future research. 
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Appendix A: Calculating the correlation length £ 

In the previous work on the model defined in (pQ), there exist both short-range and infinite- 
range interactions, but no anisotropy. Therefore we could exclude the contribution of the 
long-range correlations uniformly from the correlation function, obtaining the correlation 
length from the following relation, 

Io /2 (c(r)-c(L/2))rdr 
J L/2 {c{r)-c{L/2))dr 

Here, c(L/2) which is proportional to 7= is a good approximation for the large-r limit of 
c(r). 

In the present paper, we adopt another method because of the strong anisotropy in the 
correlation function. In a d— dimensional system, we can estimate £ from 

(k 1 -") = £ T^TT S ( k ) « e 2 '- 3 ^- (A2) 

In the present model defined by with very weak elastic interactions, there exist rather 
large, but finite clusters at the critical point due to the short-range interactions. In (1A2j) . 
we sum the structure factor with importance l/|fc| d_1 in order to avoid divergences and to 
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collect the contributions from length scales on the order of L, using the asymptotic form of 
the correlation function ( ITT)) . 



1 T. Nakada, P. A. Rikvold, T. Mori, M. Nishino, and S. Miyashita, Phys. Rev. B 84, 054433 
(2011). 

2 S. Miyashita, Y. Konishi, M. Nishino, H. Tokoro, and P. A. Rikvold, Phys. Rev. B 77, 014105 
(2008). 

3 K. L. Ronayne, H. Paulsen, A. Heifer, A. C. Dennis, J. A. Wolny, A. I. Chumakov, 
V. Schiinemann, H. Winkler, H. Spiering, A. Bousseksou, P. Giitlich, A. X. Trautwein, and 
J. J. McGarvey, Phy. Chem. Chem. Phys. 8 4685 (2006). 

4 P. Giitlich, A. Hauser, and H. Spiering, Angew. Chem., Int. Ed. Engl. 33, 2024(1994), and 
references therein. 

5 J. A. Real, H. Bolvin, A. Bousseksou, A. Dworkin, O. Kahn, F. Varret, and J. Zarembowitch, 
J. Am. Chem. Soc. 114, 4650 (1992). 

6 O. Kahn and C. J. Martinez, Science 279, 44 (1998). 

7 J. F. Letard, P. Guionneau, L. Rabardel, J. A. K. Howard, A. E. Goeta, D. Chasseau, and 
O. Kahn, Inorg. Chem. 37, 4432 (1998). 

8 A. Hauser, J. Jeftic, H. Romstedt, R. Hinek, and H. Spiering, Coord. Chem. Rev. 190-192, 471 
(1999). 

9 N. Shimamoto, S. Ohkoshi, O. Sato, and K. Hashimoto, Inorg. Chem. 41, 678 (2002). 

G. Molnar, A. Bousseksou, A. Zwick, J. J. McGarvey, Chem. Phys. Lett. 367, 593 (2003). 

1 M. Nishino, K. Boukheddaden, Y. Konishi, and S. Miyashita, Phys. Rev. Lett. 98, 247203 
(2007). 

2 C. Enachescu, L. Stoleriu, A. Stancu, and A. Hauser, Phys. Rev. Lett. 102, 257204 (2009). 

3 M. Nishino, C. Enachescu, S. Miyashita, K. Boukheddaden, and F. Varret, Phys. Rev. B 82, 
R020409 (2010). 

4 Y. Konishi, H. Tokoro, M. Nishino, and S. Miyashita, Phys. Rev. Lett. 100, 067206 (2008). 

5 I. R. McDonald, Chem. Phys. Lett. 3, 241 (1969). 

6 N. Huby, L. Guerin, E. Collet, L. Toupet, J. C. Ameline, and H. Cailleau, T. Roisnel, T. Taya- 
gaki, and K. Tanaka, Phys. Rev. B 69, 020101 (R) (2004). 

19 



K. Binder, Phys. Rev. Lett. 47, 693 (1981). 

E. Brezin and J. Zinn-Justin, Nucl. Phys. B 257, 867 (1985). 

E. Luijten and H. W. J. Blote, Int. J. Mod. Phys. C 6, 359 (1995). 

G. Kamieniarz and H. W. J. Blote, J. Phys. A 26, 201 (1993). 
X. S. Chen and V. Dohm, Phys. Rev. E 70, 056136 (2004). 

W. Selke and L. N. Shchur, J. Phys. A: Math. Gen. 38, L739 (2005). 

Y. Konishi, Graduate thesis, Graduate School of Science, The University of Tokyo, Japan (2007). 

T. T. Wu, B. M. McCoy, C. A. Tracy, and E. Barouch, Phys. Rev. B 13, 316 (1976). 

K. Boukheddaden, S. Miyashita, and M. Nishino, Phys. Rev. B 75, 094112 (2007). 

J. Stewart, O. Pohland and J. M. Gibson, Phys. Rev. B 49, 13848 (1994). 

W. Nicolazzi, S. Pillet, and C. Lecomte, Phys. Rev. B 78, 174401 (2008). 

W. Nicolazzi, S. Pillet, and C. Lecomte, Phys. Rev. B 80, 132102 (2009). 

H. Ikeda and K. Hirakawa, Solid State Commun. 14, 529 (1974). 

P. Guionneau, F. Le. Gac, S. Lakhoufi, A. Kaiba, D. Chasseau, J.-F. Latard, P. Negrier, 
D. Mondieig, J. A. K. Howard, and J.-M. Lager, J. Phys. Condens. Matter 19 326211 (2007). 
S. Lakhloufi, P. Guionneau, M. H. Lemae-Cailleau, P. Rosa, and J.-F. Latard, Phys. Rev. B 82, 
132104 (2010). 



20 



